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Abstract. We study the phenomenon of evolution by horizontal mean cur- 
vature flow in sub-Riemannian geometries. We use a stochastic approach to 
prove the existence of a generalized evolution in these spaces. In particular 
we show that the value function of suitable family of stochastic control prob- 
lems solves in the viscosity sense the level set equation for the evolution by 
horizontal mean curvature flow. 



1. Introduction. 

In Euclidean spaces, the motion by mean curvature flow of a hypersurface is a geo- 
metrical evolution such that the normal velocity at each point of the hypersurface 
is equal to mean curvature at that point. Unfortunately, even smooth surfaces can 
develop singularities in finite time, so a weak notion of evolution is necessary. 
A well-known example in is 

the surface given by two huge spheres, smoothly 
connected by a long straight cylinder. When that surface evolves by mean curva- 
ture, it splits, in finite time, in two connected components which are topological 
spheres. At the moment when this happens, the surface intersects itself. In the 
points of the self-intersection, the normal and so the mean curvature are not defined 

(e.g. m)- 

Different notions of generalized evolutions have been introduced in order to study 
the evolution of surfaces beyond the formation of singularities. Brakke defined 
in [Br| a varifold-based concept of weak solution which provides existence but no 
uniqueness, while a variational approach was developed by Almgren, Taylor and 
Wang |ATWj and Luckhaus and Sturzenhecker [LS| . Another approach, the so- 
called barrier solutions, has been introduced by E. De Giorgi |DG| and developed 
by G. Bellettini, M. Paolini and M. Novaga |BP2| IBNj . while a more recent ap- 
proach was found by G. Barles and P.E. Souganidis, [BSj . 

The notion that we are going to use follows a nonlinear-PDE-approach, found for the 
first time in 1991, independently, by Chen-Giga-Goto f |CGGj ) and Evans- Spruck 
( [ESj ). Roughly speaking, the idea consists in associating a PDE to a smooth hy- 
persurface evolving by mean curvature fiow such that the function (of space and 
time) which solves this PDE has level sets which evolve by mean curvature flow. 
Then one can define the solutions of the "generalized evolution by mean curvature 
fiow" as the hypersurface given by the zero-level sets of the viscosity solution of 
this PDE. In this paper we study the corresponding evolution in sub-Riemannian 
geometries with the help of stochastic control methods. 

Sub-Riemannian geometries are degenerate Riemannian spaces where the Riemann- 
ian inner product is defined just on a sub-bundle of the tangent bundle. To be more 
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precise, we will consider Xi, ...,Xm smooth vector fields on R" and a Riemannian 
inner product defined on the distribution 7Y generated by such vector fields. We 
assume that the vector fields verify the Hormander condition which means that the 
associated Lie algebra is equal to R" at any point. This condition has many impor- 
tant consequences. The main one is that we can always connect two points by an 
admissible path. i.e. an absolutely continuous curve such that the velocity belongs 
to the distribution for almost any time (Chow's Theorem). Therefore we can always 
define an associated distance d{x,y) on the whole R". These spaces are topologi- 
cally equivalent to the Euclidean R" but the metric is not. In fact, instead of the 
equivalence, we have only the following inequality c\x — y\ < d{x,y) < C\x — J/I"^, 
for some c, C > constants and fc > 1. For some recent applications of these 
geometries we refer to the papers by Citti and Sarti on the visual cortex ( |CS| ). 
In a sub- Riemannian space it is possible to define intrinsic derivatives of any order 
taking the derivatives along the vector fields Xi, ...,Xm- That allows us to write 
differential operators like Laplacian, infinite-Laplacian etc, using intrinsic deriva- 
tives. Wc call these operators "horizontal" . In particular we can define a notion of 
horizontal mean curvature flow. 

While there arc many results for evolution by mean curvature flow in the Euclidean 
setting, only little is known in these degenerate spaces. This evolution in a sub- 
Riemannian manifold is very different from the corresponding Euclidean motion, 
in particular because of the existence of the so-called characteristic points, which 
are points where the Euclidean normal is perpendicular to the horizontal space and 
so not admissible. Even if their geometrical meaning is very different, at the level 
of the level-set PDE, they look almost like the Euclidean "singularities," i.e. like 
the points where the Icvel-sct function has vanishing gradient. While in the Eu- 
clidean case the only "problematic" value of the gradient is zero, the corresponding 
set of gradients in our case is space-dependent and has, at each point, nonzero 
dimension. The different nature of these degeneracies creates serious difficulties 
in applying most of those techniques which are known to work for the Euclidean 
setting. To avoid the problems created from the presence of these singularities, we 
will use a stochastic approach for showing existence of solutions. 

A connection between a certain stochastic control problem and a large class of geo- 
metric evolution equations, including the (Euclidean) evolution by mean curvature 
flow, has been found by Buckdahn, Cardaliaguet and Quincampoix (in |BCQ| ) and 
Soner and Touzi (in [ST21 [ST3] ). The control, loosely speaking, constrains the in- 
crements of the stochastic process to a lower dimensional subspacc of R" , while the 
cost functional consists only of the terminal cost but involves an essential suprcmum 
over the probability space. It turns out that the value function solves the level set 
equation associated with the geometric evolution. Moreover, one can show that the 
set of points from which the initial hypersurface can be reached almost surely in a 
given time by choosing an appropriate control coincides with the set evolving by 
mean curvature flow. This stochastic approach generalizes very naturally to sub- 
Riemannian geometries. Instead of constraining the Euclidean Brownian motion, 
we use an intrinsic Brownian motion associated with the sub-Riemannian geome- 
try. This allows us to obtain certain existence results in general sub-Riemannian 
manifolds which were previously unknown. 

Our construction of controlled paths yields an analogue to the processes considered 
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for the Euclidean case in |BCQ| IST2|, IST3| , which could be called locally codimen- 
sion one constrained Brownian motion. In the Euclidean case, any control i^{s), 
taking values in the space of co-rank-one orthogonal projections induces a locally 
codimension one constrained Brownian motion B'^ as solution of the following Ito 
SDE dS" = v{s)dB. In the present sub-Riemannian case, in order to define the 
locally codimension one constrained or unconstrained Brownian motion, some ex- 
tra care has to be taken due to the geometry. We define an "horizontal Brownian 
motion" as the stochastic process whome generator is the horizontal Laplacian op- 
erator Aq = X^I^i -^i- We would like to remark that, unlike the Euclidean case, the 
horizontal Laplacian has in general also a first order part, coming from taking the 
derivatives of the vector fields. The construction of the associated unconstrained 
horizontal Brownian motion by means of the following Stratonovich SDE, is natural: 
d£,{s) = 1 ° dB\s), where B = (B\ . . . , B'") is a standard Brownian 

motion in R™. 

Replacing in previous Stratonovich SDE the unconstrained Brownian motion B by 
a locally codimension one constrained Euclidean Brownian motion B'^ in R™, we 
get the locally constrained codimension one horizontal Brownian motion £,^{s) as- 
sociated to Aq and i'(s), which constitutes a controlled path for our problem. 
Note that, in the sub-Riemannian case, the process defined by the Stratonovich 
SDE above differs from the one defined by the corresponding Ito integral just in 
the "vertical direction" (a consequence of [DGNj ). so that their horizontal projec- 
tions coincide. Moreover, in the Heisenberg group, it is possible to show, by an 
explicit calculation, that they are the same also in the "vertical" direction. 
The value function associated to this stochastic control problem is defined as 
the infimum, over the admissible controls, of the essential supremum of the fi- 
nal cost g (at some fixed terminal time T > t), for the controlled path (^'^ start- 
ing from X at the time t. We can show that this value function is a representa- 
tion formula for the generalized evolution by horizontal mean curvature flow, i.e. 
solves in the (discontinuous) viscosity sense the equation —vt + H{x, Dv, D^v) = 
where H[x, Dv, D'^v) = Aqw — Aq^oo, with terminal condition g. That means that 
u{t, x) := v{T — t, x) is a viscosity solution of ut -\- H[x, Du, D^u) = with initial 
condition g, which is exactly the level set equation of the evolution by horizontal 
mean curvature flow. So V{t) = {x € R" \ u{t,x) = 0} is a (discontinuous) gen- 
eralized evolution by horizontal mean curvature flow in general sub-Riemannian 
manifolds. Whenever comparison principles hold, we are able to show that such 
an evolution is also continuous. Unfortunately, comparison principles are still an 
open problem in sub-Riemannian geometries. As far as we know, there exists just a 
recent preprint by Capogna and Citti ( |CC| ) where the authors obtain comparison 
principles in Carnot groups but only for a particular class of initial data. In the 
same paper, the authors show also the existence, using a definition of generalized 
motion by hoizontal mean curvature flow slightly different from our notion (but 
equivalent at least in the Euclidean case). Their approach is very different from 
ours. In fact, they get the solution not by a representation formula but as limit of 
solutions of suitable approximating non-degenerate parabolic PDEs. 
The organization of the paper is the following. 

In Section 2 we introduce sub-Riemannian geometries. We recall the definitions of 
the Hormander condition and the special case of Carnot groups. The main example 
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for a Carnot group is the Heiscnberg group hence, from time to time, we will focus 
on such a particular geometry. Then we give the definition of horizontal mean cur- 
vature and of the characteristic points. Note that, there are very few interesting 
hypersurfaccs without characteristic points (for example in the Heisenberg group, 
any compact surface topologically equivalent to a sphere, has at least one charac- 
teristic point). Therefore existence results for short times with "smooth" initial 
data are not as helpful as in the Euclidean case. 

We conclude the section looking closer at the case of the Heisenberg group and 
giving several explicit examples in this geometry. 

In Section 3 we introduce a notion of generalized evolution by mean curvature flow, 
following the level set formulation introduced by Chen- Giga- Goto in [CGGj for 
the corresponding Euclidean evolution. So we associate to the evolving surface a 
degenerate parabolic equation, where the degeneracy arises where the horizontal 
gradient vanishes, as it happens at the characteristic points. At these points, the 
differential operator as a function of first and second derivatives at the point, is 
not continuous any more. In the Euclidean case this situation arises only when the 
gradient is zero, i.e. a set of dimension zero in R"^" x R". In the sub-Riemannian 
case, however, this set is both of non-zero dimension and depend on space, which 
makes the analysis far more difficult. 

In Section 4 we define and study a stochastic control problem, whose associated 
value function solves in the viscosity sense the level set equation for the evolution 
by horizontal mean curvature flow, introduced in section 3. We introduce a family 
of (Stratonovich) stochastic ODEs driven by a "horizontal constrained Browinian 
motion" and we will show that the associated generator is exactly the horizontal 
Laplacian. Moreover we study some properties of value function. In particular, 
we show a property which implies that the time evolution of its zero level set de- 
pends on the terminal value only through the zero level set, and, moreover, that 
in Carnot groups the value function is bounded and continuous in space whenever 
the terminal cost is bounded and uniformly continuous and has a a limit as |a;| — *■ oo. 

In Section 5 we show that the value function is a bounded and lower semicontinu- 
ous viscosity solution of the level set equation for the evolution by horizontal mean 
curvature flow in the sub-Riemannian case. We first derive the PDE solved by the 
value function, assuming more regularities for the solution. This proof gives also a 
justification on why the optimal control is, at any point, the projection on the hor- 
izontal tangent space of the level set. We give several examples of sub-Riemannian 
geometries covered by our existence result. We conclude dealing with the conti- 
nuity of the stochastic representation formula found. We show that if there exist 
comparison principles for the degenerate parabolic PDE introduced in section 3, 
the value function is continuous in any sub-Riemannian geometry and, therefore, 
it is a classic viscosity solution in the sense of Crandall and Lions. Nevertheless, as 
we already remarked, comparisons are known just in very few cases. 

In the Appendix we give some of the technical proofs omitted in Section 5. In 
particular we will give a rigourous viscosity proof of the main existence result. 
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2. Mean curvature in sub-Riemannian geometries. 
2.1. Sub-Riemannian geometries and Carnot groups. 

In this section we recall briefly what sub-Riemannian geometries and Carnot groups 

are (for more information, see |Be[ IHei IM] ). 

Let Xi{x), ...,Xm{x) a family of smooth vector fields on M" and 

7^, = Span(Xi(a;),...,X„(x)) 

then H ^ {{x,v) \ X G M", v £ Hx} is a distribution on M". 

Definition 2.1. A sub-Riemannian metric in M" is a Riemannian metric (•, ■) ^ 
defined on the fibers of a distribution Ti.. 

An absolutely continuous curve 7 : [0, T] ^ K" is called horizontal, if and only if, 
7(t) e H^(t), a.e. t G [0, T], i.e. 

m 

7W = E"^W^'(^W)' a.e.tG[0,T]. (1) 

1=1 

For any horizontal curve, we defined a length-functional as 

Jo 

with |7(t)|g = (7(0i7(0)|- From now to on, we choose the Riemannian metric 
{■,-)g such that the vector fields Xi, X„i are orthonormal, that means, by H]) 

1(7)-^ ^al{t) + ... + al{t)dt, 

Once defined the length-functional we can introduce the following distance 

d{x, y) := inf{/(7) | 7 horizontal curve joining x to j/}. (2) 

Let m{x) = dim(Span(7ia;)). If m{x) = n, at any point x G R", the distribution 
generates the whole tangent bundle so the Riemannian metric (t)^ induces by 
(12) a Riemannian distance on M". Otherwise, if m < n at some point, this is not 
possible and ^ can be infinite for some pairs {x,y). The so-called Hormander 
condition, which we will explain below, guarantes that d{x, y) remains finite. 
Recall that the Lie bracket between two vector fields X and Y is defined as the 
vector field which acts on smooth real functions by [X, F]/ = X{Yf) — Y{Xf). 
Let £i = {Xi,...,X,„}, £2 ^ {[X^,XJ]\i,j = l,...,m} and, for k > 2, ^ 
{[Yi, Yj] I G C'\ Yj e C\h,l^l, k - 1}\ UtTi^ C\ then the Lie algebra associ- 
ated to the distribution Ti is the set £ = IJ jgpj £* . 
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Definition 2.2. We say that the distribution 7i satisfies the Hormander condition 
if there exists r > 1 so that C = UI=i '^'^'^ span(£) = M" at any point. The 
number r is called step o/ the distribution. In such a case, we call the function 
([2]), induced by the distribution Ti, a sub-Riemannian (or Carnot-Carathcodoryj 
distance on M" , and the triple (M" , 7i, (•, is called a sub-Ricmannian geometry. 

The main consequence of the Hormander condition is that the associated sub- 
Riemannian distance is finite (Chow's Theorem). Moreover, the Hormander con- 
dition imphes that the associated distance is continuous in M" with respect to 
the EucUdean topology, and it imphes the existence of minimizing geodesies. The 
geodesies in these geometries are usuahy not locally unique. 

Carnot groups are particular sub-Ricmannian geometries, where a structure of Lie 
group is defined. We recall briefly the main definitions. For more details we refer 
to jCDPTtlDGNl lHe]. Let G = (M", •) a Lie group and g the stratified Lie algebra 
of the left-invariant vector fields. 

Definition 2.3. A Carnot group is a Lie group G = (R",-), nilpotent and simply 
connected, whose Lie algebra g admits a stratification, i.e. it can be written as 
g = (Bl^iVi, where Vi satisfy the property 

V,+i = [Vi,V,]^% fori = l,...,r-l and [Fi,K]={0} 

M" endowed with the distribution Ti. ~ Vi and the Euclidean metric on it, is a 
sub-Riemannian geometry with step r. 

On a Carnot group we can define a family of dilations 6\{x) ^ A'x, whenever 
X G Vi. Note that ||(5A(a;)||pp = A||.T||pp, where ~ d{x,0) is the Carnot- 

Carathcodory norm defined on G. Moreover, in any Carnot group, it is possible to 
define an homogenous norm, that is 



El 



with X ~ [xi, ...,Xk) and Xi € Vi. This norm is indeed homogeneous w.r.t. the 
dilations defined on the Carnot group and it is equivalent to but more useful than 
the Carnot-Caratheodory norm. In fact, explicit calculations are easy and, if we 
define the homogenous distance, setting 

do{x,y) = \\y^^ ■ x\\^ 

where is the inverse of y w.r.t. the group multiplication •, then dg is smooth 
for a = 2r!, where r is the step of the distribution. By using the dilations and the 
group operation one can get far stronger results than in the general sub-Riemannian 
case. In particular, it is possible to define intrinsic regularizations by inf- and 
sup-convolutions (see [Wj . Definition 3.1 and Proposition 3.3) which imply many 
uniqueness results for visocity solutions of nonlinear equations which are still open 
in the general case of Hormander vector fields. 
We would like to conclude the section with some examples. 

Example 2.1 (Heisenberg group). The most important sub-Riemannian geometry 
is the (1 -dimensional) Heisenberg group which is the Carnot group {E? , ■), with 
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law 



{x, y, t) ■ {x\ y', t') = (^x + x',y + y' , t + t' + 



and [x,y,z) ^ ~ {~x,~y,—z). Set p = [x,y,z) G R'^, the associated Lie al- 
gebra given by the left-invariant vector fields is X{p) = y /2)^ and Y(j)) = 
(0, 1, The bracket relations are [X, Y]=T = (0, 0, 1)^ and [X, T] = [Y, T] = 
0. Hence is a sub-Riemannian geometry with step r = 2. 
The differential of the left-translations is Lp : . q ^ p ■ q m 

/I \ 

dLp^ \ 10= -dL~^ 

The family of dilations is given by 5x{x, y, z) = (Ax, Ay, \^ z) and the homogeneous 
norm is \\ (x, z)\\^ ~ {[x^ + ?/^)^ + z^) . 

In general the n-dimensional Heisenberg group H" is defined on R^"+'^ and the 
group law, for any point (x, y, z) € R" x R" x R, is given by (x, y, z) ■ (x', y', z') = 
(x + x', y + z + + \{x ■ y' — y ■ x')) and all previous notions are still true. 

Example 2.2 (Grusin plane). The Grusin plane is the sub-Riemannian geometry 
defined on R^, by the distribution spanned by the two vector fields Xi{x, y) ~ (1, 0)"^ 
and X2{x,y) = (0,x)"^. In this case the step of the distribution is 2 (as in the 
Heisenberg group ) but the dimension of the distribution is not constant since it is 
m ~ 1 at the origin and m = 2, otherwise. Even if the Grusin plane is not a Garnot 
group, its structure is not so different from the structure ofM^. In fact, it is possible 
to define dilations and an homogeneous norm, which are Sx{x,y) = (Ax, A^y) and 
\\ix,y)\\, = \x\ + \y\i. 

Example 2.3 (Roto-translation geometry). The roto-translations geometry is gen- 
erated on R^ by Xi{x,y,9) = (cos0, sin0, 0)^ and X2{x,y,9) = (0,0,1)'^. It is a 
2-step sub-Riemannian geometry which was introduced by Gitti and Sarti in |CSj . 
in order to study the modal and amodal perceptual completion of the visual cortex. 

2.2. Horizontal mean curvature. 

Wc introduce the notion of horizontal mean curvature in sub-Riemannian mani- 
folds. There are various ways how to define mean curvature in such spaces. For 
more details and the links with the variation of the area and the approximation 
by corresponding Ricmannian objects, we refer to |CDPTj in the Heisenberg group 
and [CC, DGN, HP in a more general context. 

Given Xi, ....,Xm smooth vector fields on R", satisfying the Hormander condition, 
we indicate with M the associated sub-Riemannian manifold (R", 7Y, (•, •)g) and we 
recall that (■, ■)g is built in such a way that Xi, X„i are orthonormal, i.e. for 
V, w horizontal vectors, 

{v,w)g = {a, (3) 

where a and /3 are the coordinate-vectors of v and w w.r.t. Xi, Xm, and (•, •) is 
the usual inner product in R™. 

First we recall that the horizontal gradient of a function u : A/ ^ R is the horizontal 
vector field defined as 

Tw(x) = {Xiu)Xi{x) + ... + {X^u)X^{x) e R" 
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From now on, wc will often omit the dependency on the point x and use the 
coordinate-vector field of 

Xu= {Xiu,..,X„,uf eM'" 

Note that 

m 

1=1 

where | • | is the eculidean norm in M™. 

Before giving the main definitions we want also to point out the following notation. 
We are interested in the study of hypersurface on sub-Riemannian manifolds so 
we have to treat two different kinds of tangent spaces: the tangent space of the 
manifold and the tangent space to the hypersurface. 

Fix a point x S M. In order to avoid confusion, we call horizontal space the tangent 
space of the sub-Riemannian manifold, denoted by H^M, while the tangent space 
and horizontal tangent space are, respectively, the Euclidean tangent space of the 
hypersurface Y, <Z M and the intersection of the Euclidean tangent space with the 
horizontal space. We indicate the latter two objects by and HT^Tj. 

Definition 2.4. Let S = {u = 0} a hypersurface in M , we call horizontal normal 
of S the renormalized projection of the Euclidean normal on the horizontal space, 
which is 

Tu 

We introduce the horizontal mean curvature as the horizontal divergence of the 
horizontal normal: 

M^):=|*(^) (3) 

Unlike in the Euclidean case, the horizontal normal to a smooth hypersurface is not 
always well defined. In fact, whenever the Euclidean normal is "vertical", which 
means that its projection on the horizontal space vanishes, then rig and hence kg 
are not defined. 

Definition 2.5. Given a hypersurface S = { u = 0} C M , we call set of the char- 
acteristic points the set of the points where the Euclidean normal is perpendicular 
to the horizontal space, that is 

char{T,) = {x £ M \ H^M C T^S} = {a: € M \ HT^X = H^.M} = {x € M \\Xu\ = 0} 

(4) 

As we will see later, the existence of characteristic points make the evolution by 
horizontal mean curvature flow much different from the corresponding Euclidean 
or Riemannian evolution. 

Let u : M" ^ E be a smooth function. In the Euclidean setting it is easy to show 
that the mean curvature of E = {x G R" : u{x) = 0} is equal to the Laplacian 
minus the infinite-Laplacian, both divided by the modulus of the gradient. We 
next recall the definition of the corresponding horizontal operators and we show 
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that this fact is still true at the non characteristic points. First we recall that the 
symmetrized matrix of second derivatives is a to x to matrix defined as 

We call horizontal Laplacian and horizontal infinite- Laplacian, respectively, the 
following second order operators: 

Xu Xu 



A^u^^Xi{X,u), Ao,ooM = ( (A-^u)* 



\Xu\ \Xu\ 

Then, as in the Euclidean and Ricmannian case, it is immediate to show that 

ko{x) = \Xu\-^{l^QU - Ao,oow) (5) 

For later use, we express all the previous objects by the matrix associated to the 
sub-Riemannian geometry, the Euclidean gradient Du and the Euclidean Hessian 
D^u. So let (j{x) be the smooth mxn matrix defined as cr(a;) = [Xi{x), ....X„i(x)]"^, 
then the coordinate-vector of the horizontal gradient can be expressed as 

Xu{x) — a{x)Du{x). 

The main point is to express the symmetrized matrix of horizontal second deriva- 
tives, uing a{x). In fact, the matrix {X^uj docs not depend on just second order 
derivatives like the corresponding Euclidean one but also on first order derivatives 
coming from the derivatives of the vector fields. To be more precise, one can write 

(X'^u)* = a{x){D'^u)a'^ {x) + A{Xi, Xm, Du) (6) 

where the matrix A is a symmetric m x m matrix defined as 

A,^j{Xi,...,X^,Du) = ^{Wx,Xj +Wx,X„Du) , for i, j = 1, to (7) 

and VxiXj is the (Euclidean) derivative of the vector field Xj w.r.t. the vector 
field Xi. Hence, it is possible to rewrite previous horizontal second order operators 
as 

m 

Aou = Ti{a{x){D^u)a^{x)) + ^ {\7x,X„ Du) (8) 

1=1 

and 

Ao,.-= (H-)P^-)-^(-))fM^-fM^) (9) 
\ \a{x)Uu\ \a{x)Uu\ / 

AtY Y n ^ '^{x)Du a{x)Du \ 
MX„...,X„,,Du)j-^^,j-^^^ (10) 

This paves the way for studying the horizontal mean curvature flow by the tech- 
niques from stochastic control theory which we explain later. 

We conclude this section looking at the particular case of the Heisenberg group. It 
is known that in general S7xiXj is perpendicular to the horizontal space (see |DGN| 
for a proof in Carnot groups). An easy calculation shows that in the Heisenberg 
group S/xiXj + VxjXi = 0, for every i,j = 1,2. So 

Aou = Tr{aix)iD^u)a^{x)) and Ao,oo^^ = /(cT(x)(i?2«)a^(.T))fM^, 

^ ' \^ ' \a(x)Du\ \a[x)Du\ 
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That makes it easier to study explicit examples in the Heisenberg group. Hence, 
let S = {{x,y,z) £ M^\u{x,y, z) = O} be a surface in H-'^, then we can explicitly 
calculate all previous quantities and, in particular, the set of characteristic points 
becomes: 

char{j:) = l^{x,y,z) e \ (u^ ~ ^u^^ + ("y + |"^) = o| (H) 

Whenever E is invariant by rotation around the z-axis, i.e. u{x) = |z| — /(r), with 
r = \J x"^ + y^, we get 



^0(^,2/, ^)=±^77T^7^tW^ (12) 

(depending if z > or z < 0). In such a case the set of the characteristic points is 

c/iar(E) - {(x,2/,z) £ E | 4(/'(r))' + = o} (13) 

Obviously the only possible solutions are r = and /'(O) = 0. 

That means that the unique possible characteristic points are the flat intersection 

with the z-axis, i.e. (0, 0,±/(0)), with Neuman boundary condition /'(O) = 0. 

Using (jl2p . it is easy to calculate the horizontal mean curvature in the following 

examples. 

Example 2.4 ([CDPT]). 

(1) : Euclidean ball. // E = {(x, y, z) e M'^jx^ + ?/2 + z^ i?^}^ 

2(4 + i?2) 

Ko — I 3-. 

yjx'^ + 2/2(4 -|_ 2;2)2 

The characteristic points are (0,0,±i?). 

(2) : Koranyi ball. If = {{x,y, z) e M.^\{x^ + y^)^ + 16z^ = R^}, 



3-\/a;2 + y2 



i?2 

The characteristic points are (0,0, ±-^). 
(3): Heisenberg ball. Let E = {{x,y,z) £ M.^\d{{x,y, z), {0,0,0)) = R"^}, 
then, using the explicit formula for the Heisenberg geodesies, we have 
r — a/x^ + 2/2 = I sin(ci?/2) and z = therefore 

1 c/2 sin(ci?) — ci? cos(ci?) 

° ~ 2 sin(ci?/2) sin(ci?/2) - (ci?/2) cos(ci?/2) ' 

The characteristic points are (0,0, ±^). 

The situation is particularly easy when there are no characteristic points. 

Definition 2.6. We call regular hypersurface any hypersurface such that all 
the points are not characteristic. 

In Riemannian geometries any hypersurface is regular, while in the sub-Riemanian 
case very few hypersurfaces are. We quote the following remark due by Roberto 
Monti. 
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Lemma 2.1. Any compact surface S C H^, topologically equivalent to the 
sphere, is not regular. 

Proof. The "hairy ball theorem" from algebraic topology states that, given a vector 
field tangent to a surface in M.^, topologically equivalent to the sphere, there exists 
at least one point where the vector field vanishes. Let us consider as vector field the 
horizontal normal vector. Assuming that the surface is regular, such a vector field is 
different from zero at any point. Now we can (e.g. using the complex interpretation 
of the Heisenberg group) just rotate such a vector by ^. This new vector field is 
still not vanishing at any point but it is tangent to the surface, which contradicts 
the topological theorem. □ 

However it is possible to find some examples of (non-compact) regular surfaces 
among the rotational surfaces {\z\ = /(-\/ x"^ + y^) |, a;^ + S [a, b]} with a, 6 G R. 
Whenever a > (i.e. r ^ in whole the surface), there are no characteristic points. 
So rotational surfaces around the z-axis are regular, provided they do not intersect 
the z-axis. That remark leads to the following examples. 

Example 2.5. Regular surfaces are: 

(1) any vertical plane ax + by ^ d, 

(2) any cylinder around the z-axis, 

(3) any torus around the z-axis. 

Let us point out that non-regular surfaces are the far more interesting examples, 
because all sphere-type surfaces are not regular (Lemma 12. ip and, moreover, the 
characteristic points are what really makes this geometry so geometrically different 
from the analogous Euclidean one. 

3. Generalized evolution by horizontal mean curvature. 

In Euclidean spaces, the motion by mean curvature flow of a manifold of codi- 
mension 1 is the geometrical evolution defined by requiring the normal velocity at 
each point of the manifold. Only few results are known for mean curvature flow 
in sub-Riemannian manifolds, i.e. for the evolution obtained by replacing all the 
geometrical objects by the corresponding horizontal quantities. In these degenerate 
spaces, such a kind of evolution is very different from the corresponding Euclidean 
motion, especially because of the existence of characteristic points, i.e. points were 
motion in the (Euclidean) normal direction is not "admissible" . 

Let us define rigorously the evolution (or motion) by mean curvature flow (MCF) in 
a sub-Riemannian geometry. We give first a notion assuming that the hypersurface 
is regular (i.e. smooth without characteristic points) and then we derive a weak 
notion holding for every hypersurface. 

Definition 3.1. For t > 0, let T{t) a family of regular hypersurfaces in a sub- 
Riemannian geometry (M",7i, (•, ■)g)- l^e say that T{t) is an evolution by horizon- 
tal mean curvature flow of the hypersurface Fq if and only if the following holds: 
r(0) = Fq, and for any smooth horizontal curve x{t) : [0,T] M" such that 
x{t) £ F(<) for all t € [0,T], the "horizontal normal velocity" is equal to minus the 
horizontal mean curvature, i.e. 



vo{x{t)) := {x{t),no{x{t)))g = -ko{x{t)) 



(14) 
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where no{x{t)) and ko{x{t)) are the horizontal external normal and the horizontal 
curvature ofT{t), calculated at the point x{t) G r(t). 

Note that {x{t),no{x{t))) ^ is weh defined since x{t) is horizontal and smooth and 
we assume that T(t) is a regular hypersurface. 

([T4|) is not sufficient to describe the evolution since, like in the Euclidean case, it 
is not defined whenever the hypersurface develops singularities (which can happen 
in the Euclidean case starting from a smooth hypersurface) and it is not defined at 
the characteristic points, which are a specific feature of the sub-Ricmannian MCE. 
We introduce a weak notion of evolution by mean curvature flow, using the level 
set approach. Such a definition was given first by Chen, Giga and Goto [CGG| and, 
independently, by Evans and Spruck [ESj . It is based on the idea of defining the 
evolution of a function u{t,x) by a degenerate parabolic PDE in such a way that 
each level set {x G M" : u{t, x) = c} evolves by mean curvature as long as it is 
a smooth manifold, see e.g. [ESj . Exploiting the fact that this PDE is degenerate 
parabolic, one can define a generalized solution, called viscosity solution. 

Next we derive this degenerate PDE for regular hypersurfaces. 

Let T{t) = {u{t,x) = c}, then the horizontal normal of T{t) at x{t) is given by 

no(a;(t)) = Tu/|Tu|. Since x{t) is horizontal and smooth, 

rri 

x{t) = ^ai{t)Xi{x{t)) = cr^(x(i))a(i), for every t e [0,T]. 
1=1 

Set a(i) = (ai(i),...,a„(<))'^ and recalling that (^x{t),j^^'^ = <^a(t),^^^, 
the horizontal normal velocity can be written as 

(x(t),no(.(i))), ^ -("W'R^)^^ = \^uri<t),Du)^ 

with (•, and (•, •)^ denoting the inner product in M™ and R" respectively. 
From now to on, we can precede similarly to the Euclidean case. In fact, x{t) £ T{t) 
if and only if u{t,x{t)) ~ c. Taking the derivative in time and using (jl4p . yields 

utit,x{t)) = - {x{t),Du{t,x{t))),^ = ~\Xu\ (i(t),no(a;(i))g = \Xu\koixit)) (15) 
It remains to use ko{x(t)) = J^^i -^i ( ^\xu\ ) ' '^hich gives 

ut = Tr {{Xhi)*) - ^(A-^u)*^, = Aou - Ao^ooU. (16) 

We want to point out that equation ([TB)) is parabolic degenerate whenever Xu = 
a{x)Du = 0. We call the points where the horizontal gradient vanishes singulari- 
ties. In the Euclidean case it is known that singularities can lead to the so-called 
fattening of level sets. We say fattening occurs when the level set has no-empty 
interior, that means in particular that the gradient vanishes in an open subset, i.e. 
the co-dimension of the level set is locally zero (see }AAG[ IBP2i IGij , for more infor- 
mation). In the sub-Riemannian geometry, singularities are related to the vanishing 
of the horizontal gradient, which happens at characteristic points. 
Note that the co-dimension of the horizontal tangent space is not zero at a charac- 
teristic point. Therefore, singularities of the level set equation in sub-Riemannian 
geometries describe very different geometrical phenomena in spite of the formal 
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analogy with the Euclidean case. 

In order to introduce a generalized motion by horizontal mean curvature, we fol- 
low the definition introduced by Chen, Giga and Goto in [CGGI for the Euclidean 
evolution and by Giga in [Gi| for generic degenerate parabolic equations. 

Recall that the structure of (fT6|) is that of a degenerate parabolic equation 

ut + F{x,Du,D^u) = (17) 

with 

F{x,p,S) = -Tr{<j{x)Sa^{x)+A{x,p)) + ({a{x)Sa^{x) + A{x,p))^^,f^ 

with A{x,p) defined in (O. For sake of semplicity, set 

S = a{x)S(j^ {x) + A{x^p), 

then we can easily calculate that the upper and lower semicontinuous envelopes of 
equation (fT7|) are 

\ \(^{x)p\ \a{x)p\/ 
-tr{S)+\ma.{S), \a{x)p\^Q 



and 



-Tr(5) + (5f|^,f|^), k(.)p|^0 
\ \a{x)p\ \a{x)p\/ 



where Amax(>S') and Amin(>5') are the maximal and minimal eigenvalues of S. 

Taking S = (X'^u)* and a{x)p = Xu, we can give the following definition for the 
generalized motion by horizontal mean curvature flow. 

Definition 3.2. Let Fq = {x G M"|uo(x) = 0} hypersurface in W\ We say that 
F(t) ~ {x Cz W^\u{t, x) = 0} is generalized evolution by horizontal mean curvature 
flow if u satisfies the initial condition u{0,x) = uo(x) and it is a viscosity solution 
of (fT6|) in the sense of |Gi| , that means u is a continuous function and 

(1) for any (p G C^(R" x (0, +00)) such that u — (p has a local minimum at 
(to, xo), then 

( ipf - Aaip + Aq^oo^P > 0, at{to,xo), if Xip{to,XQ) ^0 

\(pt~ + Xma^iiX^ (fi)*) > 0, at{to,xo), if X(p{ta,XQ) 

(2) for any ip G C^(R" x (0, +00)) such that u — (p has a local maximum at 
{to, Xo), then 

( V3t - AoV3 + Ao,ooV < at{to,xo), if Xip{to,xo) ^ 

\ (ft- Aof + Amin((A'2.^)*) < 0, at{to,Xo), if Xip{to,Xo) = 



As we will see later, we can give the same definition for discontinuous functions, 
requiring the subsolution condition (Def. 13. 2[ (2)) for the upper semicontinuous 
envelope of u and, the supersolution condition (Def. 13.21 (1)) for the lower semi- 
continuous envelope of u. 
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We would like to point out that the level set approach gives a well-posed notion of 
evolution, provided that the set T{t) does not depend on the chosen parametrization 
uq but just on the level set Tq. This is the case if whenever 

U+ := {x e M" I uo{x) < 0} C {x e M" | voix) < 0} V+, 

indicating by u{t,x) and v{t,x), respectively, the viscosity solutions of equation 
P6p with initial conditions uq and vq, then 

{x e M" I u{t, x) < 0} C {x e M" I v{t, x) < 0} 

(similarly for the nonnegativc level sets). In the Euclidean case, it is possible to 
prove this, using comparison principles for the level set equation (and a suitable 
reparametrization for the initial data such that we obtain uq > vq see |Gi| , Theorem 
4.2.8. for a complete proof). 

The proof holds for equations with F strong geometric, provided a comparison 
principle for viscosity solutions holds. In our case F is strong geometric, hence 
the well-posedness of Definition 13.21 depends mainly on the existence of comparison 
principles. Unfortunately, very little is known about comparison principles (and 
hence uniqueness in the case of the evolution by horizontal mean curvature. In 
|CC| , Capogna and Citti prove comparison principles in Carnot groups for special 
classes of initial data: 

Theorem 3.1 ( [CC| . Theorems 3.1). LetG be a Carnot group andu andv hounded 
viscosity suhsolution and supersolution of ([T6|) . respectively, with initial datum uq 
and Vq. Suppose uq or vq are uniformly continuous and for any (xh^xv), (xh, Uv) ^ 
G = ffijLiVi, where we indicate by zh G Vi the horizontal part of a point and by 
zy the vertical part (i.e. zy € ffii=2^jA then uo{xh,xv) < vo{xH,yH)- Then 
u(t, x) < v{t, x), for any i > and x G G. 

The assumption on the initial datum is used to avoid the problems created by the 
characteristic points. They can circumvent the degeneracy of the equation in the 
spirit of the Euclidean proof, using the method of doubling variables and showing 
that the horizontal gradient of the chosen test function does not vanish at the max- 
imum point. Spheres, tori and any compact surfaces (see Example 12.41 for some of 
those) are not covered by the result of |CC| . 

Let us point out that the definition introduced in [CC| looks slightly different from 
ours. Indeed Capogna and Citti follow the definition of viscosity solutions for de- 
generate parabolic equations used by Evans and Spruck in [ES| . that is a more 
general notion of solution. So comparisons for the Capogna-Citti's definition imply 
comprisons for our solutions. The two definitions are equivalent in the Euclidean 
case (see jCij ) . while this equivalence is not clear in the sub-Riemannian case. We 
would like to remark that the results proved for general nonlinear degenerate par- 
abolic equations in |Gij (like equivalence of the definitions, comparison principles, 
existence, etc.) rely on techniques which are not applicable in our case. The main 
difference between the usual degenerate parabolic equations and the level set equa- 
tion for the evolution by horizontal mean curvature flow is that equation (|16p is 
discontinuous at the points (p, x) S M" x R" such that a{x)p = 0, which is a space- 
variable-depending set which has non-zero dimension in p. 

To conclude this section we are going to have a closer look at the case of the 
Heisenberg group. Since {X'^u)* ~ a{x){D'^u)a^ (x) is a symmetric 2x2 matrix. 
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h - (20) 



there are exactly two eigenvalues. As the trace is the sum of the eigenvalues, we 
can rewrite Definition 13.21 at singular points as 

- Ao^ + Xm..i{X^^)*) = (^t - A„,i„((A'V)*) > 

and 

- Ao^ + A„,in((A'V)*) - <Pt - A„,ax((A'V)*) < 

In the particular case of rotational surfaces around the z-axis in the Hciscnbcrg 
group, it is easy to see that the level set equation is continuous up to the set of 
characteristic points. Let Tq = {{x,y,z) e | |z| = f{r)}, i.e. u{t,{x,y,z)) = 
\z\ — f{t,r) with r = + ■ The level set equation at the non characteristic 
points with < r <C 1, is equal to 

A {f{r)f+ r^f"{r) 

We already remarked that the characteristic points correspond to r = with van- 
ishing Ncuman condition /'(O) = 0. So by the first-order Taylor expansion for the 
function f{r) in (that is f'{r) = /'(O) + /"(0)r + o(r) = /"(0)r), we can deduce 

lim^-/"(0) 

Hence, at the characteristic points, the matrix {X'^u)*{r — 0) has an eigenvalue 
equal to /"(O) with multiplicity 2. So in such a case the level set equation is 
continuous and we obtain ft — /"(O) = 0, whenever r = 0. In particular, there is a 
rich class of examples for which the velocity in the characteristic point is non-zero. 

4. Controlled diffusion processes. 

Let us first recall some elementary facts from stochastic analysis for continuous 
semi- martingales which can be found in any standard textbook such as e.g. |KS| . 
Given a probability space (il,JF, P) together with a filtration {J-t}t>o let £,{t) be 
continuous and adapted (i.e. ^(t) is J^t-mcasurablc), and let B{t) be a Brownian 
motion adapted to the filtration. Then the Ito integral £^dB{t)is defined as the 
following limit (as the step size of the partition decreases) in L^{^1) : 

t 2 N 

as)dB{s) := hm J2^{U){B{U+,)-B{U)). 

N^ + oo ^ — ^ 

1=1 

Note that this holds actually in a far more general setting: The convergence holds 
in the space of continuous square integrable martingales, the deterministic partition 
may be replaced by one constructed via stopping times, the integrand ^ need not 
be continuous, but merely previsible, and the Brownian motion as integrator can 
be replaced by any square-integrable continuous (semi-) martingale ri{t). In latter, 
more general case, we write {(,dr]){t) for the Ito-integral. 
The Stratonovich integral ^ o drj is defined as 

nt N 



^^+°° , = 1 



both integrals are related by the formula 
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where 77) denotes the quadratic covariation of the processes ^ and rj which is 
defined as 

Jo iV— + C50 .^^ 

The chain rule looks classical if we use the Stratonovich integral. In fact, for any 
smooth /, the process f{^{t)) satisfies 

d[fm)]=f'm)od^, 

which can be re-written as 

d[fm)] - f'mm + lf"{m)d{i,o- 

Note that, whenever ^ = _B is a Brownian motion, we have ci(^, = d{B, B) = dt 
and the formula above is the well known Ito formula. This establishes the ba- 
sic connection between second order PDE and stochastic processes which yields 
an extension of the classical method of characteristics to the case of second order 
equations. 

We would like to point out that we will use the Stratonovich calculus for defining 
our controlled stochastic processes since, because the chain rule is the classical one, 
it does not depend on the chosen parametrization and so it is intrinsic in Riemann- 
ian and sub-Riemannian geometries (see e.g. [H]). Nevertheless the Ito calculus 
will be very useful for proofs and computations (sec Sec. 5). 

4.1. The stochastic control problem. 

It is well known that viscosity solutions of certain second-order equations are closely 
related to the value function of stochastic control problems, see e.g. IFS| . The 
relation between solutions of degenerate equations like in Definition 13.21 and sto- 
chastic control problems is more complicated. Nevertheless, Soner and Touzi (in 
|ST2| [ST3] ) and, using another approach, Buckdalm, Cardaliaguet and Quincam- 
poix (in [BCQD derived a stochastic representation for a set evolving by mean curva- 
ture fiow (in the Euclidean case) . The following stochastic optimal control problem 
([24|) has very much in common with its Euclidean predecessors |ST2( IST3| |BCQ| , 
where now we have to replace the Ito by a Stratonovich SDE which reflects the fact 
that we do not work in an Euclidean space. 

Let (17, JF, {jFj}t>o, P) be a filtered probability space and i? is a m-dimensional 
Browinian motion adapted to the filtration {J-t]t>a- We define the set of admissible 
controls by 

V = {{v{s))s>Q predictable | v{s) G S*™, > 0, - > 0, Tr(/„ - j/^) ^ 1}. 

Under suitable assumptions, each i'{s) determines a (unique)control path ^*'^'''( ) 
as a solution to the SDE 

^^t,^-(-)(s) ^ V2a^(^*'^'''(-)(s)) o dB^is), s G {t, T] 

dB-'is) = L'{s)dB{s), se{t,T] (21) 

^*'^'''(-)(t) =a; 
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where odB'^ denotes the mtegral w.r.t to B'^ in the sense of Stratonovich. Using the 
relation ^ o drj = £_drj + ^(^, rj) between the Stratonovich and the Ito formulation, 
we get the following equivalent Ito formulation for SDE (PT|) 

rn 

+ J2('''(')h^^^^^(^''^'''^'^i'))ds, seit,T] (22) 

where VxiXj = DXj ■ Xi is the (Euclidean) derivative of the vector field Xj 
in direction Xi. A straightforward application of Ito's formula gives for smooth 
bounded u : R" ^ R that 

m 

duie'^'-'is)) = A/2^X,(«)(^*'-'''(s))i/(s)dB(s) 

n n 



k,l=l 



fc=l 



{e''^'''is))ds, 



where we used the notation Xi ~ {Xf, . . . , X") £ 
so that the previous identity can be written as 



Uk 



1^ and Uki 



duiC-^^-'^is)) = V2j2Mu){e''''''{s)Hs)dB{s)+tr [{,.is))^{X^u)] (e*'^'''(s))ds 

i=l 

From now on, we assume that the matrix (t{x) as well as the drift 

m 

are Lipschitz in x. Under the Lipschitz condition, classical results for stochastic 
ODEs give that for any fixed control i^, (PT|) has a unique strong (see e.g. |YZj . 
Chapter 1, Corollary 6.1). Recall that the main difference between the notions of 
strong and weak solutions is that, in the first case, the filtered probability space 
(f2, {^t}t>o, P) and the Brownian motion B are fixed while a weak solutions 
mean that there exists a process on some filtered probability space equipped with 
an adapted Brownian motion which satisfies the equation, for more details see Defi- 
nitions 6.2 and 6.5, YZJ . This difference becomes very important for the stochastic 
control problem, i.e. when considering 

infE[/(e*'-''^(T))], (23) 

where usually / is a suitably regular teminal cost function and ^*'^'''(-) are solu- 
tions of a controlled Ito SDE as e.g. (f^^ . It is clear that the properties of the 
previous minimum problem depend on the set where we take the infimum. Hence 
we define the set A of all the weak-admissible controlled pairs f |YZ| . Definition 
4.2) which are, roughly speaking, 6-tuple tt = {n,T,{Tt}t>o,^,B{-),v{-)) such 
that (^*'^'''(-), (r2,^,{:^i}t>o,P)) is a weak solution of the controlled SDE ((2T|) . 
w.r.t. the control v and the Brownian motion B in the filtered probability space 
JF, {jFt}t>o). Under certain structural assumptions for the control set and as- 
suming sufficient regularity of the coefficients, the existence of an optimal control 
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for a large class of problems as in (|23p is known, (see for example Theorem 5.3 in 
[YZ| ) . For these results it is crucial to use the weak formulation. 

Following |ST21IST3l|BCQ| , for a given bounded and uniformly continuous function 
g : R" R, we define the value function associated to the stochastic control 
probems ((2T|) . as 

Vit, x) = inf ess sup g{^'-''-'''--\T){uj)) (24) 
'^^■^ ujen 

where the set A is the set of the weak-admissible controlled pairs, defined above. 
From now to on we will often omit the dependency on uj. 

In the Euclidean case (i.e. a{x) = Id) the value function ([24|) is the solution of 
the level set equation for the evolution by mean curvature flow (backward) in the 
viscosity sense (cf. |BCQ| , Theorem 1.1), 

~Vt ^AV~ Aoo^, cc X M", t e [0, T) 

V{T,x)=g{x), xe«". ^^^^ 

Our goal is to show that this result is still true in the general sub-Riemannian case. 
We conclude this section which was devoted to stochastic control problems with a 
remark. 

Remark 4.1 (Non-Lipschitz coefficients). // the coefficients of the matrix a{x) 
and of the drift part fJ-{x) are smooth but not globally Lipschitz, the solutions of 
the SDE could explode in finite time. There are results on non-explosion for some 
classes of non-Lipschitz coefficients, but we will not investigate this issue further, 
but instead assume global in time existence of solutions for the controlled SDEs. In 
many important examples, e.g. in the Heisenberg group, the Lipschitz condition is 
satisfied and so the non-explosion follows. In particular, in this case the drift part 
is zero, so the Stratonovich and the ltd formulations coincide. 

4.2. Properties of the value function. 

In this section we study the main properties of the value function as defined by 

(IMl). 

Lemma 4.1 (Comparison Principle). Let .91,32 be bounded and uniformly contin- 
uous functions with gi < g2 on [0,T] x R", and let Vi, i € {1,2}, be defined as in 
Il24\ l with gi as terminal cost. Then 

Vi{x,t) <V2{x,t), on[0,r]xR" 
The proof is obvious and therefore omitted. 

Lemma 4.2 (Value function is geometric). Let g be bounded and uniformly con- 
tinuous, and let Vg be defined as in {2^^ with g as terminal cost. Let 1^ : R — > R be 
continuous and strictly increasing. Then 

ip{Vg{t,X)) = V^l^g){t,X) 

Proof. As Lp is increasing and continuous, (/3(inf A) = inf Lp{A) for any bounded set 
A C R. Hence, for any measurable function / : ^ R, it is trivial to note: 

(^(esssup/) = esssup(</7(/)) 

and so we can easily conclude the proof. □ 
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Remark 4.2. Lemmas \4-1\ and \4.2\ allow us to conclude, reasoning as in [Gij . that 
the sublevel set {V(t,x) < 0} depends only on the set {g{x) < 0}, and not on 
the specific form of g. Therefore we could introduce a new weak notion for the 
evolution by horizontal mean curvature by defining the zero-level set of the value 
function V{t^ x) as generalized solution at time t, which is well-defined. Moreover, 
in Theorem \5.1[ we will show that V{t, x) solves (in the viscosity sense) the level set 
equation for the evolution by horizontal mean curvature flow. Unfortunately, with- 
out comparison principles and the resulting uniqueness for solutions of the PDE, 
we cannot say that this new definition is consistent with the classical evolution for 
regular hyper surf aces. 

Lemma 4.3 (Boundedness). Assume that g is bounded and that there exists a 
weak solution of the controlled SDE (|2ip for at least one control I'is), then the 
value function V{t, x) defined in \24^ is bounded. 

Proof. The property follows immediately once we know that the infimum is taken 
over a non-empty set. □ 

In order to investigate the continuity of the value function, we have to restrict our 
attention to the case of Carnot groups. 

Lemma 4.4 (Continuity in space). Let G = (R", •) be a Carnot group, and suppose 
g : G ^ M is bounded and uniformly continuous on the one- point- compactification 
ofG, i.e. it is uniformly continuous on G and there exists lim^^^^^ g(x). 
Then V(t, x) defined in {2^^ is continuous in space. 

Proof. Denote by La{-) the left translation in the Carnot group by the element 
a € G. . As G is a Lie group, we may assume that the vector fields Xi, . . . , Xm are 
left-invariant, i.e. 

X,{a ■ x) = X,{La{x)) = {DLa){X,{x)), (26) 

for i = 1,...,TO, where DLa is the derivative of the left translation (see e.g. [W] 
for more details on Lie groups). Let ^*'^'''( ) be a constrained codimension one 
horizontal Brownian motion, with d^*'^''^^'^ = \/2a'^ {^*'^''^^'\s)) o dB^{s), then, by 
the chain rule for Stratonovich integrals, it holds 

d [la (e*^^^"-' ))) = (DLa) o (V2a^ (e*'^'^"(s)) o clB'^is)) 

^V2{DLa) (a^ (C*'"'"^ ^*))) ° dB^is) = (i, (e^'^'^^-^s)) ) o dB'' {s) 

(27) 

where we have used (|26p for the last equality. Hence the left translation of a 
codimension 1 horizontal Brownian motion yields another one. Now fix a point x, 
e > and choose a control v.j. such that 

V{t,x)^t> esssupg{C''''''^^\T)). 

Let a = y ■ x~^. By (P7|) . the path ) starting at the time t in y, is equal 

to Lai^*''^''^''^'^). (Note that the control i^x is the same for both points x and y.) 
Therefore 

V(t,y) < esssupg(77*'^'^-(-)(T)) = esssupg(La(C*''''''^^-Hr))) 

= esssup(.g(^*'^'''-(-)(T))-h (g(ia(C*''^'''^^'HT'))) 

< V{t,x) + e + esssu-plg {LaiC''''''''^-\T))) ~ g '-'^ {T))\ 
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Choose a large number R > then 



ess sup 



9 (?*'^'''-"(r)(c.)) - g [Ue''^''^^^-HT)H)) 



< sup \g{z) - g{a ■ z)\ + sup \g{z) ~ g{a ■ z)\ =: A + B 

{zeR": \z\<R.} {zGR": \z\>R} 

where we set z = ^*'^'''-(-)(T)(tj) and so a • z = La{^^''''''-'^-\T){uj)). 
Note that |a-x| — > oo if|a;| —^ oo. Therefore we can use the continuity of g at 
oo to find a sufficiently large R such that B < e. As, by continuity of the group 
operation, |a • cc — a;| ^0 (uniformly on compact sets) as \a\ = \y ■ x^^\ 0, we 
can use the uniform continuity of g to find S > such that V{t, y) < V{t, x) + 3e 
for |x — j/| < 5. Reversing the role of x and y yields the continuity. □ 

5. Existence of a generalized evolution by horizontal mean 
curvature flow in sub-rlemannian manifolds. 

Using the value function for the stochastic control problem introduced in the pre- 
vious section as representation for the viscosity solution of equation p6|) . we get 
an existence result for the generalized evolution by horizontal mean curvature flow 
as given in Definition 13.21 By classical results (see e.g. }FS| . [T]), it is known how 
to find the equation solved by value functions of the form inf,yg^E[5(^*'^''^'')(T)] . 
Unfortunately, the value function V{t,x) defined in ([24| looks different, because of 
the essential supremum instead of the expectation. Hence the idea (already used 
in |BCQ| ) is to approximate formula ((24)) with functions that look like the infimum 
of an expectation and then to pass to the limit, essentially using the fact that the 
ip-norm of a fixed nonnegative function converges to the essential supremum as 
p oo. 



Since we are only able to show that the value functions defined in ([24|) are lower 
semicontinous, we need to use the viscosity theory for discontinuous functions. Next 
we will recall the definition. For more details on this theory, we refer to (Baj . 

Definition 5.1. A locally bounded function u : M" x [0,r] ^ R is a discontinuous 
viscosity solution of equation p5|. ifu*{t,x) is a viscosity suhsolution and Uf{t,x) 
is a viscosity supersolution of the same equation, where u* and are respectively 
the upper and lower semicontinuous envelope of u, i.e. 

u*{t, x) :~ inf{u(t, x)\v cont. and v > u} ~ limsup{ii(s, y)\ \y ~ x\ < r, \t ^ s\ < r}, 

uJt^ x) := sup{u(i, x)\v cont. and v < u\ = liminf{it(s, y)\ \y — x\ < r\t — s\ < r}. 

The main result of this paper is the following existence theorem. 

Theorem 5.1. Let g : M" M bounded and Holder continuous, T > and 
(j{x) = [Xi{x), Xm,{x)]'^ amx n-Hormander matrix with m < n and smooth 
coefficients. Assuming that cr(x) and '}2^iLi'^ XiXj{x) are Lipschitz (in order to 
have non-explosion for the solution of the SDE), then the value function V{t,x) 
defined by (|24p is a bounded lower semicontinuous viscosity solution of the level set 
equation for the evolution by horizontal mean curvature flow (P5)) . with terminal 
condition V{T,x) = g{x). 
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Remark 5.1 (Examples). We recall that the Lipschitz requirement of the theorem 
is always satisfied if the vector fields Xi, Xm and their derivatives are smooth 
with at most linear growth a infinity. Hence the previous theorem covers many sub- 
Riemannian geometries, as, for example, the Heisenberg group, the Grusin plane 
and the roto-translation geometry, introduced in Examvles \2.1[ \2.S\ and, \2.3l 

Next we give a proof that works rigourously just under additional, strong assump- 
tions on the regularity. In fact, for passing to the limit, we will need a strong 
convergence of the approximating functions together with all the derivatives in- 
volved. In general those assumptions are not satisfied and indeed V{t, x) is just 
a viscosity solution and not a classic solution. In the Appendix, we will present 
a viscosity proof that follows the one used in the Euclidean case by Buckdahn, 
CardaHaguet and Quincampoix in |BCQ| . Nevertheless, the following proof gives a 
clear idea of why the optimal control, realizing the infimum in (|24p . is given by the 
projection on the tangent space of the level set of the solution at that point. 
Before we start with the heuristic arguments, let us introduce a regularization which 
will be used both for the heuristics and for a rigorous viscosity proof. As the es- 
sential supremum is the limit of L^'-norms, we consider, following |BCQ| , for any 
1 < p < +00, the following value functions: 

Vpit,x)=M^ (E[.9P(e*'^^''(■)(^))])'^^ (28) 

Note that we can replace g hy ag + b for real numbers a and b and therefore, as g 
is bounded, assume C > g{x) > 1, (for any C > 0) 

The idea is to derive the PDE solved by the value functions ([28|) and then to show 
that V is their limit as p +00 and solves a limit-equation which is exactly the 
level set equation for the evolution by horizontal mean curvature flow, (|16[) . In fact, 
we have 

Lemma 5.1. Under the assumptions of Theorem \5.U we have 



V{t,x) = lim Vp{t,x) for all {t,x) g [0, T] x {pointwise convergence) 

p — v+oc 

(29) 



with V(t,x) as in {24^ . 

As the norms are increasing and bounded by the essential supremum, for each 
fixed control, it is clear that V{t,x) > \imVp{t,x). In order to show equality, we 
can argue as in |BCQ| , noting that in our case the controlled SDE has a drift part, 
but it depends on the control only through i/^, and our control set is, as the one in 
|BCQ| , convex in i/^. Hence similar arguments, based on |YZ| . apply. Let us now 
explain heuristically how to find the equation solved by V. 

Heuristic proof. We first look at 

Up{t,x) = V^{t,x) = inf E[5P(e*'^'^(T))] (30) 

It is known (see e.g. jFSj . [T] ) that Up{t, x) is a viscosity solution of 
{ -{Up)t+H{x,DUp,D^Up) = Q, a;GR", te [0,T) 
1 Up{T,x)^gP{x), xeR" 



where 



H{x,p, S) = sup 



^h:v[F{x,y)F^{x,y)S]+{pi{x,v),p) 



(31) 
(32) 
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whenever is a solution of the stochastic ODE 



In order to compute the structure of the Hamiltonian associated to ([30]) , we have to 
consider the Ito formulation of the stochastic control problem, given in (P^ . that 
means F{x,v) = ^/2a'^ {x)v{s) and ^,{x,v) = Y^^j^ii^'^i^))^ x^Xj{x). Hence ((5^ 
becomes 



H{x,p,S) 



sup 



-TT{a{x)Sa'^{x)v'^{s)) 



{v\s)\^,{Wx.X,{x),p) 



(33) 



for any x,p Cz M" and any symmetric n x n matrix S. 



Unlike for Euclidean mean curvature flow, our Hamiltonian depends on the gradi- 
ent and on space, not just on S, (In fact in the Euclidean case Vxi^j = and 
(j{x) — Id). Assuming g locally Lipschitz (or locally Holder) and cr(x) smooth 
(hence locally Lipschitz), it is possible to show the continuity of Up, exactly as in 
Propositions 1.3 and 1.4 in [T]. So Up is a (continuous) viscosity solution of (PT|) 
with H given by (|33p . 

Once we know the equation solved by Up{t,x) = V^{t,x), we can easily find the 
corresponding PDE for Vp. 

Since C > g{x) > 1, > 1 > 0, too, and we can divide by pV^~^ . Assuming that 
all functions involved are smooth, a trivial calculation tells that Vp solves 



f -{Vp)t + H{x, DVp, {p - l)Vp-'DVpiDVpf + D^Vp) = 0, a; € M", i € [0, T) 
\ VpiT,x) ^ g{x), xeW 

. . 

Whenever Vp is just continuous, we can show that it solves equation p4p in the 
viscosity sense, by applying the previous calculation to the (smooth) test functions. 
The continuity for Vp follows from the continuity for Up. 

The following computations conclude the proof, assuming Vp and V are so regular 
that 



{Vp)^ Vt, DVp DV, D^Vp D^V as p ^ +oo 



As we already remarked, those assumptions are in general not satisfied and so we 
will give in the Appendix a more rigourous viscosity proof. 

However, the following computation explains why the supremum in (j33|) is attained 
whenever V ^ Im — no (g) no, where no is the horizontal normal to the level set (at 
the non-characteristc points). 

So, let us first assume XV{t,x) — a{x)DV{t,x) ^ (imphes XVp{t,x) ^ 0, at 
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least, for large p). We can write explicitely the Hamiltonian in ([34|) : 
Hix, DVp, (p - l)V„^DVp{DVpf + D^Vp) = 



sup 



-{p- l)Tr[V-^vv'^{a{x)DVp)[a[x)DVpf] + Tv[uv'^ (j{x)D'^Vpa'^ {x)] 



(35) 



Recalling that (X^Vp)* ~ cr{x)D^Vpa^ (x) + A{x, DVp) where A{x,p) is defined by 
([7]), we observe that 

Tr(A'Vp)* = TT{a{x)D^Vpa^{x)) + Tt{A{x, DVp)) 

The trace of the first order part is 

m 

TrA{x, DVp) = J2 {^x.X,{x), DVp) 

i=l 

so Hamiltonian (|35p can be rewritten in "horizontal notation" as 
Hix, DVp, (p - l)Vp-^DVp{DVp)'^ + D^Vp) = 



sup 



-{p- l)Vp^TY[vy'^ {XVp){XVp)'^] + 1:y[vv'^ {X'^Vp)*] 



(36) 



Note that 



~{p-l)Vp-^TT[uv^{XVp){XVp)^] ^ -{p-l)Vp-^TT[{v^XVp){v'^XVp)^] < 0, 

and so it goes to —oo as p — *■ +00. Hence, in order to attain the supremum, we 
need (at least for large p) that XVp = 0. Since the horizontal gradient is in the 
direction of the horizontal normal, the optimal control 77 has to coincide with the 
projection on the tangent space, that means 



V = Irn -no®no 



(37) 



To get the level set equation, we have to write Hamiltonian ([55]) in a bit different 



way. Let I„i be the m x m identity-matrix, we can replace v by /„ 
a e M™, then, for any n x n matrix S, it holds 



with 



sup I 



TT{iy^S)] 



sup I 



Tr((/™ ~a(E) a)S)] = -Tt[S] + max {Sa, a) 



Using the optimal control ([37|) and recalling that no = XV/\XV\ and S = (X'^V)*, 
we can conclude that the limit Hamiltonian, as p — > 00, is 



H{x,DV,D^V) = -Ti[{X^V)*] 



, , XV XV 
(X^V)* 



\xv\' \xv\ 



So the limit equation of problems (|34[) . as p ^ +00, is exactly the level set equation 
for the horizontal evolution by mean curvature, at non-characteristic points. 
Next we would like to investigate what happens at the characteristic points, so 
we assume XV{t,x) ~ 0. Passing to the limit in the first order disapears 
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whatever the control v looks like. Then we get, for any control v = Im — a® a and 



= -Ft + sup 



= -Vt ~ Tr(A'^F)* + maxTrfaa' {X'-V)*) 

|a| — 1 

= -Vt - AoF + A,„ax(A'V)* 
So we find, as expected, the upper semicontinous regularization of the equation. 

We like to conclude this proof with some remarks on the regularity of x) and 
the initial condition. 

First we can note that V(t^x) is bounded since the datum g is. Moreover, since 
is a non-decreasing sequence of continuous functions, we have 

V(t,x)= lim K(t,a;) = supK(t,a;), (38) 

Hence Vit^x) is, a priori, just lower semicontinuous. To conclude, we observe that 
F(r, x) = limp^+oo 14'(T, x) = g{x). □ 

Example 5.1. In the case of the Heisenherg group, all the hypersurfaces quoted 
in Examples \2.4\ are covered by the existence result, given by Theorem \ 5.1l In 
fact, whenever the hypersurface is compact and (like in the Euclidean ball), we 
can choose as initial datum g, the signed Euclidean distance in a neighborhood, 
extended continuously by constants outside, so that the corresponding g is bounded 
and Lipschitz continuous. For the Carnot-Caratheodory ball, we can define g{x) = 
d{x, 0) — 1 locally around the surface and constant outside. The datum in this case 
is bounded and locally Holder continuous (since it is Lipschitz w.r.t. the Carnot- 
Caratheodory metric d). The same holds for the Kordnyi ball. 



We would like to show that the functions defined in ([24|l are continuous in time 
and space and to get so the existence of a continuous viscosity solution of the level 
set equation for horizotal mean curvature flow. 

Unfortunately, so far we are able to show directly the continuity of the value function 
((24)1 just w.r.t. X and for Carnot groups (Lemma |4.4[) . 

Following the strategy introduced in |BCQ| , it is possible to get the continuity of 
the solution via comparison principles for viscosity subsolutions and supersolutions. 
The problem with generalizing this idea to the sub-Riemannian setting is that 
comparison principles are known just for special initial data in Carnot groups, 
hence so far we are able to get a full continuous existence result just for the class 
of initial data covered by Theorem 13.11 by [CC| . 



Proposition 5.1. Let g : E" — > R bounded and Holder continuous, T > and 
cf{x) an mx n-Hormander matrix like in Theorem \5.1\ If comparison principles for 
(|25p hold, then the value function V{t,x) defined by (|24p is the unique continuous 
viscosity solution of the level set equation (j25p . satisfying V{T,x) = g{x). 

Let V'^{t, x) the half-relaxed upper-limit of Vp{t, x), i.e. 

V^{t,x)= limsup Vp{s,y) (39) 
{s,y) (t,x) 

p +00 
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We remark that, in our case V'^{t, x) = V*{t,x), (more details on this point, will 
be given in the Appendix). We can in this way prove the following lemma, which 
is the key-point in order to apply comparison principles and obtain the continuity 
of the solution. 

Lemma 5.2. For any x e M", V'^{T,x) < g{x). 
For a proof of the lemma, we refer to the Appendix. 

Proof of Proposition \5.1[ By Theorem 15.11 and Lemma 15.21 once we have compar- 
ison principles, it is trivial to conclude. In the viscosity proof of 15.11 (given in the 
Appendix), we will show that V*{t,x) = V'^{t,x) is a viscosity subsolution while 
V^{t^x) = V(t,x) is a viscosity supcrsolution of equation (|25p . Since, by Lemma 
15.21 V^T,x) < g{x) while g{x) = V{T,x), comparison principles imply V'^{t,x) < 
V{t,x). Moreover V'^{t,x) > V{t,x) by definition. Hence V'^{t,x) = V{t,x), which 
means V{t, x) upper semicontinuous. Since V{t, x) is already lower semicontinuous 
as supremum of continuous functions, wc conclude that V(t,x) is continuous. □ 

The previous conditional result becomes a full result in the cases where comparison 
principles are known, as in some particular hypersurfaces in Carnot groups ( |CCj ). 

Corollary 5.1. In Carnot groups and under the assumptions of Theorem \5.1\ and 
Theorem \3.1\ the value function V{t^x) defined in ([24]) is the unique continuous 
viscosity solution of the level set equation (|25p. with V{T,x) = g{x). 

6. Appendix. 

In this appendix we give the proofs omitted in previous section. 

Next proof uses ideas from |BCQ| , applied to the horizontal case, to give a formal 

viscosity proof of the main theorem. 

Proof of Theorem \5.1\ We have to show that V{t, x) solves the horizontal level set 
equation in the sense of Definition 15.11 

First we recall that, since V{t,x) is lower semicontinuous in time and space, then 
V^(t,x) = V{t,x). Let V^{t,x) the half-relaxed upper-limit of Vp(t,x), defined in 
(p9| . We remark that V'^ > V and upper semicountinuous. Since the upper semi- 
continuous envelope V* is the smallest upper semicontinuous function above V, then 
V*{t,x) < V^{t,x). Moreover since Vp{t,x) is non-decreasing, Vp{t,x) < V{t,x), 
for any x, t and p > I. Taking the limsup in x, t, p, we get also the reverse inequality 
V^{t,x) < V*{t,x), hence V* = VK Therefore, to verify Definition ISH we have to 
show that V(t, x) is a viscosity supersolution and V'^{t, x) is a viscosity subsolution. 

First we show that V{t, x) is a supersolution of So let ip G C^([0, T]; C2(R")) 
such that V ~ if has a local maximum at (t, x). Two different cases occur. 
If Xip{t, x) ^ 0, we have to verify that 

- iftit, x) - AoV3(t, x) + Ao,oo¥'(i, > (40) 

while, if X(p{t, x) = 0, we need to check 

-(^t(t,a;)- Ao¥>(t,x)-|-An,ax((A'V)*(i,a;)) >0 (41) 
Note that, for any p > 1, there exists [tp, Xp) such that Vp — Lp has a local minimum 
at {tp,Xp) and (tp,Xp) — > {t,x). In fact, we can always assume that {t,x) is a 
strict minimum in some Bii(t,x). Set K = Bn,(t,x), the sequence of minimum 
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points {tp,Xp) converge to some £ K. As 1^ is the limit of the Vp and lower 

semicontinuous, therefore a standard argument yields that {x,t) is a minimum, 
hence it equals {x,t). Since is a solution of then 

-iptitp,Xp) + H{xp, {p-l)Vp~^Dip{Dipf + D^ip){tp,xp) > 

In the case a{x)D(p{t, x) ^ 0, we can write the Hamiltonian in a more explicit way. 
Set 

Si^ {p- l)Vp^ Xip{tp, Xp)iXip(tp, Xp))"^ 
S2 = {X^(p)*{tp,Xp) 

H{xp, 81,82) = -Tr(S'i + S'2) + Amax(S'i + '5'2) = -Tr(S'i)-Tr(52) + Amax(S'i + S'2) 
= -{p - l)Vp-\tp, Xp)\Xip{tp, xp)\^ - Ao^itp, Xp) + Xn,USi + 82) (42) 

since the trace operator is linear and Ti-{{X(p{xp)){X(p)'^ (xp)) ~ \X(p{xp)\'^. 
Now we need the following result. 

Lemma 6.1 ( |BCQ| , Lemma 1.2). Let 8 a symmetric m x m-matrix such that the 
space of the eigenvectors associated to the maximum eigevalue is of dimension one. 
Then, 8 Amax(5') is in a neighborhood of 8. Moreover, Z?Amax('S')(iJ) = 
(^Ha,a'), for any a £ R™ eigenvector associated to Amax(»S') and \a\ ~ 1. 

We apply previous lemma to the matrix 

^ Vp'^{t,x)Xip{tp,Xp){Xip{tp,Xp))'^ 

which clearly satisfies the assumptions of previous lemma. 

Expanding the Hamiltonian (j42p around 8p and then, passing to the limit as p — > 
+00, we get exactly ([iO]). 

The case XLp{t, x) = is much easier. We have just to use the subadditivity of the 
function 8 — > Amax('5') and remark that, since Vp is supersolution 

0<-ipt+ H{xp, DVp, {p - l)Vp-^Dip{D^f + D^ip) 
< -^t -{p- l)Vp-'\X^f - Tr((A'V)*) + A„,ax((p - l)Vp-'X^{X^f + (X^y) 

< ^ipt -{p- 1)^-11^-^12 _ ^ _ l)Vp-^\Xip\^ + A,„ax((-%'V)*) 
at the point [tp,Xp). So, passing to the limit as p —t +00, we find (|41[) . 

To verify the subsolution property for y* = y«, let tp £ C^iO, T]; C2(R")) such 
that V'^ ~ (f has a maximum at (tp, xo) ^^'^ we may assume that such a maximum 
is strict. Let {tp,Xp) a sequence of maximum points oi Vp — ip, we can find a 
subsequence converging to x. Hence, since Vp are solutions of (j34[) . we have 

< -<^t + H{x, {p - \)ip-^ Dip{Dipf + DV) (43) 

at {tp, Xp). We define, for any z > 0, x, d £ M" and any n x n symmetric matrix 8. 

Hp{x,z,d,8) = - ~ '^^ \a[x)d + A{x, d)\^ - Tria^ {x)8<j{x) + A{x, d)) 



z 



^ — ^-^{a{x)d){a{x)df + a^{x)8a{x) + A{x, d)] (44) 
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and 

! -Ti{<j'^ {x)Sa{x) + A{x, d))+ 
({a^{x)S<j{x) + A{x,d))-^^^,.^^^, \<j{x)d\ ^0 
- Tiia^ix)Sa{x)) + A,nax(a^(a;)5a(x)), \a{x)d\ =0 

(45) 

Since Amax(o'"^5CT + A) > Amin(o'^S'(T + A), it is clear that 

7?, {x, d, S) = -ff* (x, d, S) 

Moreover, as in |BCQ| , we can observe that 

Hp(x, z, d, S) > H{x, d, S), for any z 

which is trivial for |a(a;)d| = (by (|44p ) and for |(T(a::)d| ^ 0, it follows taking 
a = i^li^j^i in the definition of maximum eigenvalue as Ainax(»5') = max|a|=i {Sa, a). 

Set z = ip~^{tp,Xp) > 0, d ~ Dif{tp,Xp), S = D^ip{tp,Xp), by (|43|) taking the 
limsup as j3 ^ +oo, we can deduce that 

> ipt+H4x,Dip,D^ip) 

at {t,x). That concludes the proof since it shows that the value function V{t,x) 
satisfies Definition 15.11 □ 

We conclude by giving a proof of Lemma 15.21 the key point for obtaining the 
continuity of the functions defined by the stochastic representation formula via the 
comparison principle. 

Proof of Lemma \5.'2\ Assuming that it is not true, there exists a point xq such that 

vHT,xo)>gixo)+e 
for e > sufficiently small. Then we use as test function 

(p{t, x) = a{T -t)+ f3\x - xof + gixo) + - 

with a > — C/3, with C a constant depending just on the data of the problem and 
the point xq (in the Euclidean case C = — 2(n — 1)) and f3 > 1 sufficiently large. 
Now we can find a sequence (tx, Xk) ~* (T, a;o) and pk —> +oo as A: — > +oo such that 
Vp^. — if has a positive local maximum at some point (sk,yk), for any fc > 1 (see 
|BCQ[ for more details). To get the contradiction we need to use the fact that Vp is 
a solution (so in particular a subsolution) of equation ([M]) and, by the subsolution 
condition, get 

a + C/3 < 
which contradicts the choice a > —CjS. 

Unfortunately in our case, unlike in the Euclidean case, the test function, inserted 
in the equation for Vp, does not give a constant number since the Hamiltonian 
depends on the space- variable. Nevertheless, we can observe that the functions Vp 
are bounded uniformly in p so, by the growth of ja; — xq], the maximum points are 
such that Uk & Bfl(cco) =: K, with R independent of k. 

Using this remark and the continuity of the terms which we get calculating H in ip 
and at the point (s^, yk), we will be able to show the following lower bound 

a + Cf3<-ipt + H{{p - l)ip-^Dip{D(pf + D^ip) < 
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and so conclude. 

First we need to remark that 

(pt{t, x) = —a, Dip{t, x) = 2P\x — xo\, D'^ip(t, x) = Id 

Remarking that at the point (s^., j/^), we have 

0>{p- l)^-^D^{D^f + D\) >a- Tr(a(yfc)a^(y,) + A{x, x - .To)) 

+ Amin(a'(x)(T^(a;) + {A{x,x - xq)) 

Recahing that there is a compact set K such that yk E K for all fc, we get by 
continuity 

> (p - l)(f-^Dip{D(ff + D^ip) >a + 2CP 

with 

C = — max Tr(CT(a;)(j"''(a::))— max A(a::, a; — a;o)+minAmin(cr(a;)o-"^(a;))+min Xniin{A{x,x—xo)) 

With such an estimate, we are able to obtain the same contradiction as in the 
Euclidean case, choosing a > —0(3. □ 
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